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Abstract 

We report a full (3D) micromagnetic simulation of a set of 100 ferrite (Fe-304) 
cylindrical dots, arranged in a 10 by 10 square (planar) array of side 3.27 /xm, 
excited by an external in-plane magnetic field. The resulting power spectrum 
of magnetic excitations and the dynamical magnetization field at the result- 
ing resonance modes were investigated. The absorption spectrum deviates 
considerably from that of a single particle reference simulation, presenting a 
mode-shifting and splitting effect. We found an inversion symmetry through 
the center of the array, in the sense that each particle and its inversion 
counterpart share approximately the same magnetization mode behavior. 
Magnonic designs aiming at synchronous or coherent tunings of spin-wave 
excitations at given spatially separated points within a regular square array 
may benefit from the new effects here described. 
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1. Introduction 



For several decades, there has been a great interest in the study of the 
collective spin excitations in magnetically ordered media, and recently the- 
oretical and experimental investigations have been thoroughly conducted 



|5]. These investigations established the prospect of controlling 
SWs in magnonic crystals (similarly to the control of light in photonic crystals 
jsj]), motivating a whole new field, currently being referred to as magnonics 
(of. 7| and references therein). 

In magnonic crystals, dipolar (magnetostatic) interactions have an impor- 
tant physical role, since they couple excitation modes of individual, closely- 



spaced partic 



net iz at ion 



es, affecting both the static and dynamic behavior of the mag- 



in the form of Bloch waves 



10|. This effect results in the formation of collective modes 

ed and forbidden 
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12 



magnonic states at given frequencies, or band gaps 15 



16] . These and other 



particular characteristics stimulated new research directions in the study of 
"spin- waves" (hereon SWs), given the possibility of designing filters and 
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waveguides for microwave nanotechnology applications 

However, only recently experimental and numerical investigations on the 
modification of normal modes of magnetic excitations in peri odically ar- 
ranged nanomagnets have been carried out in some detail 
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26|, |27|, given the need for advanced experimental and computational 
capabilities. The general theoretical formulation of magnetic phenomena at 
scales of ~ 10~ 6 — 10~ 9 m is based on the Landau-Lifshitz-Gilbert (LLG) 



equation 
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30] for the magnetization dynamics. Note that the LLG 



equation only be solved analytically for special cases 31], hence computa- 



2 



tional micromagnetics with increasingly detailed simulations have been an 



important aid at understanding microma gne t ic p henomena 



32 



33]. 



In a series of papers by Kruglyak et al. |23|, [24|, |25|, |26| , particular attention 
was given to the investigation of the magnetization dynamics of square arrays 
of submicron elements of different sizes under a range of bias fields. These 
investigations involved the use of time-resolved scanning Kerr microscopy 
to probe the magnetic response of nanoelements, along with micromagnetic 
simulations to aid the analysis of the resulting spectra. These works have 
generally shown that the position of the mode frequencies as a function of the 
element size as well as their relative absorption amplitudes present a com- 
plex dependence and follows a nonmonotonic behavior. It was also observed 
that the position of mode frequencies did not follow the prediction of the 
macrospin model for an isolated element. It was inferred that a nonuniform 
distribution of the demagnetizing field could be responsible for nonuniform 
precession within the elements of the array, adding to the complex depen- 
dence the role of exchange interaction. It has also been noted that, as the size 
of the element in the array is decreased, the edge regions of a given element 
present increasingly dominant modes confined by the demagnetizing field in 
relation to uniform modes. In addition, in studies where the orientation of 
the external magnetic field was rotated in the element plane of the array, 
the effective magnetic field inside a given element also presented an "extrin- 
sic" anisotropic contribution due to the stray field from nearby elements, 
as contrasted to an "intrinsic" anisotropy occurring in an isolated element. 
Furthermore, a dynamical configurational anisotropy was necessary to quali- 
tatively explain the data. An additional important feature, specially observed 
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23| 



in micromagnetic simulations of arrays of nanoelements, is the splitting of 
precessional modes, a feature experimentally verified recently as detectable 
collective spin wave modes extended throughout the array 

Here we report a 100-particles micromagnetic simulation - to our knowl- 
edge, the largest detailed micromagnetic calculation ever performed at the 
given scale and number of particles, with a careful observation of the accuracy 



requirements (see accuracy details in Section 2.3 of Ref. 34 L which were also 
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34j motivated 



adopted here). This work is part of an ongoing project 
by the theoretical investigation and design of new nanostructured magnetic 
configurations with interesting SWs or magnonic band gap behavior, suitable 
for different applications in the microwave frequency range. We show that 
the collective magnetization behavior is constrained by an inversion sym- 
metry through the center of the array. In particular, this opens interesting 
possibilities for applications that require spatially coordinated patterns. 

2. Methods 

The simulations were performed by using the freely available integrator 



OOMMF (Object Oriented Micromagnetic Framework) 36], which was in- 
stalled and executed on a 3 GHz Intel Pentium Desktop PC, running Kuru- 
min Linux. The present 100-particles simulation turned out to be a computa- 
tionally demanding one (taking about 4 months for completion, considering 
interruptions), and no other variation of the parameters were attempted at 
ihis time. We applied the same methodology described in our previous works 



35 



34j, based in the procedure given by Jung et al. [37] 



An incident in-plane magnetic field was applied uniformly to the ferrite 
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particles, composed by a static (dc) magnetic field (B dc = fioH dc ) of 100 
mT in the y direction, and a varying (ac) magnetic field (B ac = fioH ac ) 
of small amplitude (1 mT) in the x direction, with the functional form: 
B ac = (1 — e ~ Xt )B ac fi cos(ut). We varied the frequency (/ = u/(2ir)) from 2 
to 9.8 GHz, in steps of 0.2 GHz, resulting in 40 different OOMMF frequency 
runs. The time domain of the applied B ac field was discretized in intervals 
of 0.005 ns. 

OOMMF performed the numerical integration of the LLG equation lead- 
ing to the evolution of the magnetization field of 100 ferrite particles regu- 
larly distributed in a square, 10 by 10 array of side 3.27 /mi. The particles 
were identical cylindrical dots (each with a diameter d = 0.3 /im and thick- 
ness 5 = 85 nm). We adopted a small interparticle (edge-to-edge) spacing 
(s = 0.03 /im <C d). The simulation was stopped at 5 ns, giving 1000 out- 
puts for each frequency run. A reference simulation of a single-particle with 
a diameter d = 0.3 /im was also conducted, with the same global parameters 
of the 100-particles simulation. Table [1] lists the parameters used to set up 
the OOMMF integrator in both cases. 

The power spectra of magnetic excitations were obtained by excluding 
data from the first 2 ns of the averaged magnetization vector in the x di- 
rection, (M) x (t < 2 ns), and taking the Fourier transform of the remaining 
time domain data, (M) x (2 < t < 5 ns). The maximum Fourier peak at 
each frequency provided the magnitude of absorption at the given frequency. 
A spline fit to the absorption data was performed in order to facilitate the 
comparison of the overall behavior of the curves, but individual data points 
were maintained. 
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3. Results 

The power spectra of magnetic excitations for the 100-particles simulation 
and single-particle simulation are shown in Fig. [TJ The resonance peak in 
the power spectrum of the single-particle simulation is clearly split into four 
distinct peaks of lower amplitude in the 100-particles simulation, with two 
peaks at a lower and the other two at a higher frequency with respect to the 
reference fundamental mode, which is located at the gap between resonances 
2 and 3 of the 100-particles simulation. 

We analysed the nature of the modes of interest by an inspection of the 
time-dependent magnetization vector field configuration. Bitmaps or "snap- 
shots" of the corresponding simulations were generated, selected at two points 
of the ac field cycle, namely, at the highest (r) and lowest (r + tt) representa- 
tive amplitudes of the equilibrium magnetization oscillation. We subsampled 
the x-component of the magnetization field in order to show 9 vectors per cell 
element, and different pixel tonalities correspond to different values of M x . 
In Fig. [TJ (inset), the magnetization vector field of the one-particle simulation 
around its resonance peak is shown. This should be contrasted with those of 
Fig. [2] (snapshots of the 100-particles simulation at the previously identified 
four peaks of interest). 

Fig. H2a) shows some zoom-out regions of the 100-particles simulation in 
order to allow for the identification of several types of magnetization mode 
behavior, which will be discussed in more detail in the next section. On the 
other hand, from a visual inspection of the 10 by 10 snapshots, it is possible 
to identify an inversion symmetry through the center of the array in such 
a way that each particle and its inversion counterpart share approximately 
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the same mode behavior in the array. In other words, by setting a Cartesian 
grid on the array, where the origin of the coordinates is the center of the 
array, and where each element is centered at coordinates an inversion 

transformation — > leaves the magnetization configuration of 

the array approximately invariant. Some examples are shown in Fig. H^b). 
This global symmetry should arise from dipolar couplings, but the exact for- 
mulation is yet not clearly understood. This effect has already been pointed 
out in our previous work in the cases of 2 by 2 and 3 by 3 array simulations 



(c.f. Fig 8 of Ref. 



35] and discussions in Ref. |34j), but it was not possible 



at that time to extrapolate whether the effect would persist in a larger array. 

In order to address in a more quantitative way the magnetization dy- 
namics distribution in the array, we computed two simple estimators, which 
nevertheless establish the relevance of the visually noted effect. The first 
estimator is the modulus of the difference of average magnetization values 
at the points of the ac field cycle: =| (M x ) ^(t) — (M x )^(t + 7r) |; 

where (M x )^ is the average magnetization of a given particle at grid coor- 
dinates The second estimator, <J(ij), is the modulus of the difference of 
standard deviation values of particle magnetizations within the array, that 
is, the standard deviation is computed with reference to the average magne- 
tization of the whole array. Fig. @] shows the resulting m^ t j) and maps 
computed for the 4 modes of interest. Notice that each map pixel is labelled 
by the particle coordinate (i, j). We discuss these maps in more detail in the 
next section. 
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4. Discussion and Conclusion 

It is presently understood that particular features of the power spec- 
tra of magnetic excitations can be associated with the various types of 



nodal 

HQ 



behavior of the time-dependent magnetization field (see, e.g., Refs. 



20 



23 



24 



25 



26{|). A general "spin- wave" behavior (SWB) for the 
excitations indicates that the magnetization vectors present small ampli- 
tude oscillations about a nonuniform static magnetization field. In addition, 



the following subclasses of excitations can be identified [20|, |34| : (i) "Quasi- 
uniform" behavior (QUB): the movement of each magnetization vector is 
similar to that of its neighbors, except for the regions around the edge of the 
particle; and (ii) "Edge-like" behavior (EDB): the magnetization field in the 
central region of the particle is almost entirely static and aligned with the 
direction of the external dc field; the nonuniformly distributed magnetization 
vectors present small amplitude oscillations near the edges of the particle. In 
particular, these modes may be more affected by the dipolar couplings from 
nearby particles. 

In the present work, the nature of the reference absorption peak is basi- 
cally due to QUB, whereas the 100-particles simulation presents all types of 
behavior (as can be seen from an inspection of the zoom-out regions exem- 
plified in Fig. [3]^a)). In particular, peak 1 at 3.8 GHz is dominated by QUB, 
with a few EDB cases specially for particles at the borders of the array. Peak 
2 at 4.8 GHz shows a mixture of QUB and EDB, with a few SWB in some 
particles at the top and bottom of the array. The number of SWB cases 
appears to increase in peak 3 at 6.0 GHz and is the dominant absorption 
mechanism of peak 4 at 7.2 GHz, specially for particles located towards the 
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center of the array. This effect is approximately the same as that observed 
previously for 3 by 3 arrays (c.f. Fig. 9 of Ref. |34j). 

In order to address in a more quantitative way the inversion symmetry 
in the array (c.f. Fig. [2]), as pointed out in the previous section, we refer to 
the mu j) and <ruj) estimators shown in Fig. HI It is clear that w,uj\ gives 
larger values for larger amplitudes of the average particle magnetization, 
but also the original inversion symmetry through the center of the array of 
the magnetization configuration should be translated approximately into a 
symmetry around a central horizontal axis in a m^j) map. It is clear that if 
the standard deviation distribution is similar at the two opposite amplitude 
points, then a map of <ruj) should be uniform. This estimator also gives larger 
values for particles whose behavior with respect to the whole array presents 
a substantial difference at the two points of the cycle, thus furnishing an 
overall measure of the degree of cycle "coherence" . 

The expected central horizontal axis symmetry is indeed seen in Fig. H] 
(left column of panels), which corroborates our visual analysis. In peak 1, 
m(jj) is larger for particles localized in the top and bottom rows of the array. 
In peak 2, note the interesting regular, alternating magnetization excitations 
along specific rows of the array (also inferred from a visual inspection of the 
snapshots). Peaks 3 and 4 show higher values of muj) in the central regions, 
in contrast to peak 1. The cruj) estimator results (Fig. HI right column of 
panels) show that, excluding the top and bottom rows of the arrays in cases 
1, 3 and 4, the array behavior is reasonably similar at the two extremes of the 
cycle. However, case 2 shows again a pattern, where the second and ninth 
rows (related to particles with very low amplitudes of the average particle 
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magnetization, c.f. left column of panels) have standard deviations that 
significantly differ at the two extremes of the cycle. 

Qualitatively, the results reported in the present work, specially the split- 
ting of the resonance mode and its decreased relative amplitude, as well as 
the spatially nonuniform behavior of the elements in the array, are compati- 



ble wit 
works 



i the behavior of the collective excitations reported in similar previous 



23 



24 



25 



. Particularly, in Ref . |27| , for the first time the col- 
lective spin wave modes of the entire array has been experimentally detected. 
With the aid of a micromagnetic simulation of a 3 x 3 array, the authors ob- 
served how the row/column elements behaved out of phase with the rest of 
the elements of the array, qualitatively explaining the nature of the splitting 
of the precessional modes of the absorption spectrum. Our results show that 
this behavior is even more complex when considering a more extended array 
of 10 x 10 elements. Yet, as already mentioned, a symmetry pattern can be 
identified. Indeed, the inversion symmetry here noted is compatible with the 
experimentally observed "tilt" of the modes in regions of higher amplitude 



27] , relatively to the horizontal and vertical axes, reported in that work. We 



believe that such a "tilt" would be observed in a larger simulation for that 
material. 

In summary, our present results and previous indications allow us to infer 
that, for interparticle (edge-to-edge) separations at least of the order of ~ 10 
to 20% of the particle's diameter, dipolar couplings in periodically arranged 
cylindrical nanomagnets will cause a global, coherent magnetization behavior 
across a square array distribution, with an inversion symmetry through the 
center of the array (for in-plane magnetic excitations). The power spectrum 
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shows a clear four-fold resonance feature, wherein the magnetization field dis- 
tribution and dynamics show interesting patterns and trends. We hope that 
the effect here reported will stimulate the development of a theory that will 
generally describe and predict similar mode symmetries in periodic arrays. 
Our results may be of interest for magnonic device architectures, specially 
for applications dealing with pairwise SWs excitations of submicron ferrite 
particles at spatially separated points in a periodic square array. 
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Tabic 1: Main OOMMF parameters, fixed for all simulations. 



Simulation Parameter/Option 



Saturation magnetization [A/m] 
Exchange stiffness [J/m] 
Anisotropy constant [J/m 3 ] 
Anisotropy Type 
First Anisotropy Direction (x,y,z) 
Second Anisotropy Direction (x,y,z) 
Damping constant 
Gyromagnetic ratio [m/(A.s)] 
Particle thickness [nm] 
Cell size [nm] 

Demagnetization algorithm type 



Parameter 



5.0 x 10 5 

1.2 x lir 11 

-1.10 x 10 4 

cubic 

(111) 

(10 0) 

0.005 

2.21 x 10 5 

85.0 

5.0 

const, in each cell 
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Fig. 1: Top panel: Power spectrum of magnetic excitations of the 100- 
particles simulation. Bottom panel: Comparison of the former spectrum 
(diamond symbols) with that of the single-particle reference simulation (dot 
symbols). Four distinct peaks in the spectrum of the 100-particles simulation 
are identified. "Snapshots" of the magnetization vector field of the one- 
particle simulation around the resonance peak are indicated by the arrow. 
Snapshot to the left refers to the highest amplitude of the oscillation and the 
snapshot to the right, to the lowest one. 

Fig. 2: "Snapshots" of the magnetization vector field of the 100-particles 
(10 by 10 particle array) simulation around the four peaks identified in the 
previous figure. For each peak, the upper snapshot shows the magnetization 
state of the array at the highest amplitude of the oscillation (r) whereas the 
snapshot immediately below, at the lowest amplitude (r + n). Symmetry 
axes are shown schematically at the top left of the figure. 

Fig. 3: (a) Amplified particular selections (3 by 4 sub-arrays given by 
the solid rectangular shown inside the 10 by 10 array snapshots), at the two 
points of the cycle, (b) Illustration of the symmetric mode behavior: one 
example of particle pairs is taken from each of the four peaks at r. Grid 
coordinates are indicated above the particles. 

Fig. 4'- Estimator maps (m^j) and a^j) distributions) for the arrays, as 
explained in the main text. New symmetry axis is indicated. 
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